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We study cosmological solutions to the low-energy effective action of heterotic string theory in- 
cluding possible leading order a 1 corrections and a potential for the dilaton. We consider the 
possibility that including such stringy corrections can resolve the initial cosmological singularity. 
Since the exact form of these corrections is not known the higher-derivative terms are constructed 
so that they vanish when the metric is de Sitter spacetime. The constructed terms are compatible 
with known restrictions from scattering amplitude and string worldsheet beta-function calculations. 
Analytic and numerical techniques are used to construct a singularity-free cosmological solution. At 
late times and low-curvatures the metric is asymptotically Minkowski and the dilaton is frozen. In 
the high-curvature regime the universe enters a de Sitter phase. 
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I. INTRODUCTION 

Arguably, the most perplexing problem of modern cos- 
mology is the initial singularity problem. The theorems 
of Hawking and Penrose state that many of the manifolds 
of General Relativity are geodesically incomplete Q] . In 
particular, many solutions to the Friedmann-Robertson- 
Walker (FRW) isotropic and homogeneous model of the 
universe contain singularities. In order to gain a com- 
plete description of the early universe a theory of quan- 
tum gravity is required. In general, it is believed that 
such a theory may somehow resolve the initial singular- 
ity allowing us to obtain well-defined, finite solutions to 
calculations of physical quantities. 

Superstring theory is currently the best candidate for 
a theory of quantum gravity. It is therefore only natu- 
ral to try to use string theory to probe the very early 
universe. Many of us arc hopeful that string theory (or 
M-theory) will lead to a consistent cosmological model 
capable of resolving the initial singularity [45|]. We are 
inspired, in part, by the ability of stringy physics to re- 
solve certain sing ularities in time-independent orbifolds 
and conifolds |10|-[T^. A full understanding of the initial 
singularity problem will most likely require nonperturba- 
tive physics, and most certainly will require a better un- 
derstanding of time-dependent solutions in string theory. 

In this paper we provide a toy model of a nonsingular 
FRW cosmology, based on string theory. Our goal is to 
discuss only one possible way in which string theory may 
begin to address the initial singularity problem. Our con- 
crete starting point is the low-energy effective action of 
heterotic st ring theory including possible leading order a' 
corrections j4g. While the exact structure of these cor- 
rections is not known, we provide an example in which 
the resulting cosmology is singularity- free (47J ■ To study 
the system we use both analytic and numerical tech- 
niques. The metric is asymptotically Minkowski space- 



time at low-curvatures and evolves to de Sitter space in 
the high-curvature regime. The significance of the a' 
corrections is controlled, in part, by the dilaton field 
These higher-derivative terms become important in 
the high-curvature regime and are constructed to allow 
de Sitter solutions in the early universe. In this way 
our model provides the first natural realization of the 
Limiting Curvature Construction (LCC) in terms of a 
well-motivated, physical theory [24|-[3(|. Our analysis 
involves a dynamical dilaton and a novel, string-inspired 
form for the higher-derivative terms. 



II. THE ACTION 

In Z)-dimensions, the string tree-level effective action 
for the massless boson sector (in the string-frame) is: 
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where $ is the dilaton, a' = £ 2 and the tilde indi- 
cates that we are using the string- frame metric c/mn, 
where M, N = 0, . . . , d. In the above, £o takes on the 
values 1/4, 1/8, for the bosonic, heterotic and type 
II superstring theories, respectively. The "• • •" refer 
to other higher-derivative terms of order a' 2 , and the 
Lagrangian C2 represents the leading-order corrections 
to the action and is made up of four-derivative terms 
C 2 {R 2 , R^W", R^\aR^ Xa i ■ ■ ■) HIIH. We choose to 
work within heterotic sting theory. While the two-point 
action for the heterotic string includes a gauge vector 
boson A n (with field strength tensor F mn ) and an anti- 
symmetric tensor H a bc- 
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we consider only the most relevant massless modes, the 
dilaton $ and graviton g mn . In order to obtain results 
that are easy to compare with General Relativistic theory 
we work in the Einstein-frame obtained via a conformal 
transformation of the metric g a p = e^g a p. Applying 
this transformation to 1|2.1|1 gives the Einstein-frame ac- 
tion 
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where £ = A/(D - 2) and A = 2(D - 10)/3a' is a con- 
tribution to the cosmological constant that vanishes in 
ten-dimensions. In the Einstein-frame the dilaton cou- 
ples directly to the higher-derivative terms e - ' 5 */^- In 
this way it is possible to control the "strength" of the £2 
terms, in part, via the dilaton. 

Our knowledge of the exact form of £2 is incom- 
plete. Requiring the action to reproduce string theory 
S-matrix elements can determine only some of the coef- 
ficients of potential covariant terms in £2- This is be- 
cause some terms do not contribute to the S-matrix or 
provide contributions that overlap in form with those 
of other terms [H, Some sort of off-shell super- 

string calculation is required to fix the exact coefficients 
of these terms. Scattering amplitude or string worldsheet 
/3-function calculations predict only the Riemann squared 
term R^\ a R^ Xa in £ 2 . In Pt Metsaev and Tseytlin 
fix the contribution 
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where L2 is the term being multiplied by a' in l|2.3ll . 
However, terms such as R 2 and R^R^ may be added 
and subtracted from C2 by performing field redefinitions 
of C-2pt H3- Because of this imprecise knowledge of £2, 
£2 is commonly assumed to be the Gauss-Bonnet invari- 
ant 
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An advantage of choosing this particular structure for £2 
is that the resulting equations of motion are second-order. 
The cost of choosing an invariant of a form other than 
()2.5|) is that our equations of motion will be fourth-order. 
Because we are interested in finding a nonsingular theory 
we will desi gn something other than the Gauss-Bonnet 
form for £ 2 [48| | . A simple way to ensure that our theory 
has a nonsingular solution is to choose an £2 that will 



vanish for a nonsingular spacetime. We can then look for 
solutions that approach this nonsingular manifold in the 
large curvature regime. As predicted by the string the- 
ory calculation (|2.4|) . we take the leading term in £2 to 
be Riiv\aR^ vXa ■ An elementary set of nonsingular space- 
times are the set of maximally symmetric spacetimes of 
constant curvature. The metrics of constant curvature 
are locally characterized by the condition (in D = 4) 
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The space with constant curvature and R = is 
Minkowski spacetime. The space for R > is de Sitter 
spacetime, which has topology R 1 x S 3 . The space with 
R < is anti-de Sitter spacetime, which has topology 
S 1 x R 3 . 

Using equation (|2.6|) it is easy to construct the desired 
invariant that will vanish for spacetimes of constant cur- 
vature (in D = 4) and whose leading term is R^uXcrR^ 
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Note that it is possible that there are other solutions 
satisfying £2 = (see [47]). This is a limitation of the 
construction as presented. If there are other solutions 
satisfying this condition then we could try to modify this 
Lagrangian so that these solution no longer obey this 
condition. However, we are not primarily interested in 
finding all the solutions to this theory. Our main goal 
is to show only that nonsingular solutions of the type 
described above exist. The four-dimensional action is 
given by equation Ij2.3|l . and may be written 
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where £'($) = a'e~ /8. In the above we have assumed 
a potential for the dilaton field V r ($) and we have ignored 
the contribution to the cosmological constant term, Ae 2 * 
[i.e. we do not attempt to solve the cosmological con- 
stant problem). We will comment further on the dilaton 
potential below. 

Variation of the action with respect to the metric ten- 
sor yields: 
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while variation with respect to the field $ gives 




(2.11 ! 



where a(i) is the scale factor of the universe and n(t) 
is the lapse function. For simplicity we take the metric 
to be flat (fc — 0). We also choose a time-dependent, 
homogeneous dilaton <E>(i). An alternative derivation of 
the equations of motion 1(2. lOj) and l|2.11|l is achieved by 
substituting the metric l|3.12(l into (|2.9|) and varying the 
action with respect to a, n, and <I>. The resulting EOM 
are: 



III. COSMOLOGICAL SOLUTIONS 



We now study cosmological solutions to the theory 
(|2.9|) . We assume a time-dependent FRW homogeneous 
and isotropic metric 
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18Aa 4 + e 2 *a 4 (t/($) + 2$ 2 ) = 6a(2Ad 2 (2d4> + a) + a(e 2 *d 2 — Ad 2 + 2Ad(— 2$d + a))) (3.14) 

12A(d 2 - ad) 2 + e 2 *a 3 (12d$ + a(V"($) + 4$)) = 0, (3.15) 



where A = a'/8 and the prime in the last equation de- 
notes differentiation with respect to <&. In the above 
we have set n = 1. Note that the dynamics of this 
system are completely determined by the Goo equation 
(|3~14l) and the EOM for the scalar field QZXffy . Equation 
(13.13(1 serves as a constraint equation. To reduce these 
equations to a second-order system, introduce H = a/a. 
Equations 13.14fl and l|3.15|l reduce to 

H 2 - 2HH - 6H 2 H + 4HH& 

e 2 * / • \ 

= -— (V^) + 2<P 2 -6\H 2 ) (3.16) 
6A V / 

12XH 2 = -e 2 * f^M+4$ + 12W) ,(3.17) 



A. The late-time universe 

At late times (low curvatures) our action must 
mimic General Relativity and the Einstein-Hilbert action 
Seh = (1/2k 2 ) J d A x^—g R . This implies the contribu- 



tion L>($)£ 2 - 2(V$) 2 - in equation (E3 must 

not produce any physical deviations from Seh m the 
low-curvature regime. Furthermore, the dilaton must 
approach a <3? i— > $ . If the dilaton remained time- 
dependent it would produce observational consequences 
such as variations in gauge-couplings. We also demand 
V($) i— > to avoid gaining a contribution to the cosmo- 
logical constant in the late universe. A de Sitter space so- 
lution a(t) — ao exp Ht corresponds to constant H = Hq. 
In the case under consideration (V($) i— > 0, $ i— > <& ) the 
EOM are satisfied by a de Sitter solution with the con- 
straint H = and hence, a(t) = const. Therefore, in the 
(late-time) case of constant dilaton and vanishing po- 
tential the equations of motion are solved by Minkowski 
spacetime (the special case of de Sitter with H = 0) with 
flat metric g a p — rj a p — diag(—l, 1, 1, 1). 
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FIG. 1: The time evolution of the potential V($). For this 
plot (and all subsequent plots) we take t\ = 4 and Ho ~ 1. 



B. The early universe 

In order to remove the initial Big-Bang curvature sin- 
gularity present in Einstein gravity at t = 0, the term £2 
must become significant. This will force the metric solu- 
tion to de Sitter spacetime with metric (|3.12[1 and a(t) oc 
exp Ht. Discovering a solution which evolves smoothly 
from a de Sitter phase at early times to Minkowski at 
late times will provide an example of a universe which is 
everywhere nonsingular. We now produce such a solution 
and proceed to study its stability against classical per- 
turbations. Simple forms for $ and V(Q>) that provide 
the desired behavior and obey the EOM are 

3(t) = $ tanh(^-^ , (3.18) 

and 

V($) = F (($ + l) 2 - 4) , (3.19) 

where we set $0 = 1, Vq = —3Hq/2 and t = 1. 
The constant t\ is taken to be large enough so that 
$(i = 0) w -1. Note that the potential (|3~Tg|) is every- 
where positive for all allowed values of the dilaton $ (over 
the range, (—1,1)) and hence, for all time t € (—00,00) 
(see Fig.l). The choices for $ and V are certainly not 
unique. Any $ and V that have the correct asymp- 
totic behavior have the potential to generate nonsingu- 
lar solutions of the type under consideration. Inserting 
the <!> ansatz l|3.18|l and the V ansatz (|3.19|) into the 
EOM (equations (|3~T3|) to (|?TH)|) ) we find that as t h-> 0, 
a(t) oc exp Ht (in our specific example Hq = 1). At late 
times the metric approaches flat Minkowski space. 



C. Numerical analysis 

While we have shown that our solution has the desired 
asymptotic behavior, namely, de Sitter space at early 



times and Minkowski space at late times, we must prove 
that this solution remains nonsingular for all values of t 
in between. To do this we conduct a numerical study of 
the equations of motion. Consider the phase space gen- 
erated by the equations l|3.14[l and (|3.15ll with $ and V 
given by (|3~T%j) and The ($,if)-plane is plotted 

in Fig. 2. Note that the range of (—1, 1) is mapped to 
the entire history of the universe t S (—00, 00). 

Clearly, the curvature is bounded for all time and 
therefore, the universe is free of curvature singularities. 
At early times ($ 1— > —1) the curvature \H\ is approxi- 
mately a positive constant corresponding to the de Sit- 
ter phase. At late times the dilaton approaches a con- 
stant («& 1— > 1), the curvature vanishes and the spacetime 
approaches flat Minkowski space. Of course, recent ob- 
servational evidence from supernovae data |4Cj and the 
WMAP satellite ^lj indicate that the universe is not only 
expanding but undergoing an accelerated expansion. It 
is easy to modify the Minkowski condition in this model 
in order to accommodate various late time behaviors (for 
example, by adding an appropriate matter Lagrangian 
to the system). To achieve an accelerating universe at 
late times one simply could tune the potential to not go 
exactly to zero at late times like in quintessence models. 

It is important to discuss the stability of this solution. 
We now show that this solution is not (in general) an 
attractor at early times but is at late times. Therefore, 
nonsingular solutions of the type described here are not 
generic in the model under consideration. It is possible, 
however, that other nonsingular solutions exist in this 
theory that are attractors. Furthermore, as pointed out 
by Starobinsky pjj, the evolution of the Universe need 
not follow a generic solution, it may well be described 
just by this unique one, at least initially. The late time 
behavior of Minkowski spacetime with frozen dilaton is 
classically stable against small perturbations. Let us be- 
gin with the stability analysis of the early-time solution. 
For our considerations it is sufficient to consider stabil- 
ity against homogeneous fluctuations. An analysis of in- 
homogeneous fluctuations is considerably more involved. 
Let 

a(t) = a e m (1 + S(t)) , \S(t)\ « 1 . (3.20) 

$(*) - $0 tanh (1 + *(*)) , \<t)\ « 1 • 

. (3 - 21) 

Recall, that when 1 1— » 0, $ 1— > — 1. To first approximation 
in S and e the Goo equation and the equation of motion 
for $ become 

X 

6 + 3H5 + — = (3.22) 
4e + 12He = 3H 2 e. (3.23) 
The solutions to these equations are 

V9Ae 2 ff 2 -4 t 

5= dl+d2e 7 Tl (3-24) 

exp(((3ffe + ^ppH) )i/2 e) 
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FIG. 2: The (<!>, iif )-phase portrait generated from the equa- 
tions of motion. Note that the curvature remains finite for 
the entire history of the universe. 



order a' corrections to the low-energy effective action. 
A nonsingular solution is given in which the universe 
evolves from an early-time de Sitter phase to a late-time 
Minkowski spacctime with constant dilaton. One limita- 
tion of this model is that at high-encrgics and large cur- 
vatures E 2 a' ~ 0(1) higher-order curvature corrections 
(e.g. a' 2 corrections) become important. Including these 
terms or quantum loop corrections in g B — could, 
presumably, re- introduce the singularity. Of course, if the 
string scale is TeV then such corrections will become im- 
portant sooner then if the string scale is near the Planck 
scale. It is likely that a full nonperturbativc analysis is 
required in order to fully understand the initial singular- 
ity problem. Our analysis is meant to show a possible 
way in which string theory may address this issue. We 
have shown that including higher-derivative terms in the 
action can (under certain circumstances) resolve the ini- 
tial curvature singularity. 



ci + c 2 e 



e ((3+2V5) Ht)/2 



(3.25) 



Here e?i, g?2, c\ and c-i are constants. The conditions 
that e and S be small (<< 1) at the beginning of our 
analysis when t ~ implies that \c\ + C2I << 1 and 
Mi + (h\ « 1- While S remains small for all time (for 
small d\ and cfe) e, in general, grows greater than unity 
and our solution becomes unstable [4!|. Hence, at early 
times the solution is stable against perturbations of the 
metric but not against perturbations of the dilaton $. 
This demonstrates that nonsingular solutions of the type 
given above are not generic. 

To see if the late-time Minkowski space solution is sta- 
ble we consider perturbations around the Minkowski so- 
lution 



Finally, it would be interesting to see if this construc- 
tion can regulate the bounce that occurs in the Pre- 
Big-Bang (PBB) model of string cosmology 0, |44| . 
The analysis would involve ideas from this paper and 
from [32I l33| . In the PBB model the universe starts 
out in the perturbative string Minkowski vacuum. The 
universe then goes through a high-curvature, collapsing 
regime during which the a' corrections become impor- 
tant. Hopefully, after a successful "graceful exit", the 
universe enters an expanding low-curvature radiation (or 
quintessence) dominated phase. In this paper we have 
studied only the post-collapse branch of the PBB model. 
The collapsing branch would be the time reverse of this. 
We leave a more detailed analysis of the PBB scenario in 
this context to a future paper. 



o(t) = a (1 + 5(t)) , |<5(i)|«l 



(3.26) 



and perturbations of the dilaton of the form (|3.21l) . Re- 
call that at late times (as t 1— ► 00), <j> 1— > 1. To first 
approximation in S and e the Goo an( i the equation of 
motion for <i> give 5 = e = 0. Therefore, our late-time 
Minkowski space solution is an attractor. 
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